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ABSTRACT 

Context. This is the second of a series of papers presenting VLBI observations of the 293 Caltech-Jodrell Bank Flat-Spectrum 
(hereafter CJF) sources and their analysis. 

Aims. To obtain a consistent motion dataset large enough to allow the systematic properties of the population to be studied. 
Methods. We present the detailed kinematic analysis of the complete flux-density limited CJF survey. We computed 2-D kinematic 
models based on the optimal model-fitting parameters of multi-epoch VLBA observations. This allows us to calculate not only radial, 
but also orthogonal motions, and thus to study curvature and acceleration. Statistical tests of the motions measured and their reliability 
have been performed. A correlation analysis between the derived apparent motions, luminosities, spectral indices, and core dominance 
and the resulting consequences is described. 

Results. With at least one velocity in each of 237 sources, this sample is much larger than any available before and allows a meaningful 
statistical investigation of apparent motions and any possible correlations with other parameters in AGN jets. The main results to 
emerge are as follows: - In general motions are not consistent with a single uniform velocity applicable to all components along a 
jet. - We find a slight trend towards a positive outward acceleration and also adduce some evidence for greater acceleration in the 
inner-most regions. - We find a lack of fast components at physical distances less than a few pc from the reference feature. - Only 
~4% of the components from galaxies and <2% of those from quasars undergo large bends i.e. within 15° of ±90°. - The distribution 
of radial velocities shows a broad distribution of velocities (apparent velocities up to 30 c). 15% of the best sampled jet components 
exhibit "low" velocities which may need to be explained in a different manner to the "fast" motions. - Some negative or "backwards" 
superluminal motions are seen and in 15 cases (6%) these are definitely significant. - We find a strong correlation between the 5 GHz 
luminosity and the apparent velocity. - The CJF galaxies, on average, show slower apparent jet component velocities than the quasars. 
- The mean velocity in the VLBA 2cm survey (Kellermann et al. 2004) is substantially higher than in the CJF survey - the ratio could 
be roughly a factor of 1.5-2. This supports the observed trend of increasing apparent velocity with increasing observing frequency. 
Conclusions. This AGN survey provides the basis for any statistical analysis of jet and jet-component properties. 

Key words. Techniques: interferometric - Surveys - Galaxies: active - Radio continuum: Galaxies 



1. Introduction 

Extragalactic radio jet sources have been intensively studied 
for several decades now, for a number of important reasons. 
These include their association with supermassive accreting 
black holes, which are recognized to play a substantial role in the 
formation and evolution of galaxies. Radio jets can help to reveal 
their energetics, and the way they interact with their environ- 
ment. The formation, acceleration, and propagation of radio jets 
also challenges our understanding of the laws of physics under 
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extreme circumstances, such as (ultra-)relativistic speeds, high 
magnetic field strengths, and, possibly, the occurrence on galac- 
tic scales of Poynting flux jets, and/or, of non-baryonic matter 
In addition, radio sources can be bright enough to be visible at 
very large redshifts, and thus they allow population studies to be 
carried out over a large fraction of cosmic time. 

The parsec-scale properties of jets have been observed with 
all available VLBI networks at all available wavelengths (e.g., 
the review of Zensus 1997, and references therein). Immediately 
alongside morphological studies have come studies of the mo- 
tions of features seen within these jets. Those motions often ap- 
pear to be faster-than-light, or superluminal. This phenomenon 
occurs as a result of light travel time compression when features 
move with a speed close to c at a small angle to the line-of- 
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sight; a useful introduction to the phenomenon can be found in 
Pearson & Zensus (1987). Detailed studies and intensive mon- 
itoring campaigns of individual objects have proven to be im- 
portant to understand jet component motions in a handful of 
selected sources (e.g., 3C120: Gomez et al. 2001; 1803H-784: 
Britzen et al. 2005a&b; 0735H-178: Agudo et al. 2006). It has 
emerged that different sources can reveal quite different motions; 
see for example the compilations of Barthel et al. (1988), Zensus 
& Pearson (1990), GhiselHni et al. (1993), Vermeulen & Cohen 
(1994), Jorstad et al. (2001), Kellermann et al. (2004), Cohen et 
al. (2007), and Piner et al. (2007). This indicates a clear need to 
study large, well-defined samples that can be subjected to statis- 
tical analysis. 

In the 1980s and 1990s it was thought that a careful study 
of the statistics of superluminal motions as a function of red- 
shift could have a major impact on cosmology (e.g., Cohen et 
al. 1988; Vermeulen & Cohen 1994) and that this would be the 
major driver for studies of large samples of sources. The atten- 
tion has now turned to the study of the astrophysics of the jets 
themselves. 

The CJF survey (Taylor et al. 1996a) on which we report here 
provides a large, well-defined, complete flux-density limited flat- 
spectrum sample that has been treated fully homogeneously as 
regards observational strategy, data reduction, and data analysis. 
Images at 5 GHz have been obtained for nearly all objects from 
global VLBI or VLBA data at least at three observing epochs 
from 1990. A detailed description of the observations and data 
is presented in Britzen et al. 2007a (hereafter Paper I). With at 
least one velocity in each of 237 sources, this sample is much 
larger than any available heretofore. We have not just derived a 
single speed per object, but have been able to obtain motions for 
several individual features in many of the sources. We have also 
studied curvature and accelerations, rather than obtaining only 
total speeds, or radial motions. 

Preliminary results from some subsets of the CJF data have 
already been published in Vermeulen (1995), Britzen et al. 
(1999, 2001b), and Britzen (2002). Polarization information of a 
sub-sample of CJF-sources has been published in Pollack et al. 
(2003). Detailed information on the VLBI observations, data re- 
duction, analysis, imaging, model fitting, and component identi- 
fication is given in Paper I; the database is available online at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/AH-A/ In addition, 
the data can be downloaded from a CJF-archive at the MPIfR in 
Bonn (http://www.mpifr-bonn.mpg.de/staff/sbritzen/cjf.html). 
The present, second paper, deals with the derivation of motions 
from the component positions. The motion statistics are ana- 
lyzed in depth, and correlations with some other parameters, 
in particular the observed luminosity, are discussed. In a third 
paper (Britzen et al. 2007b, Paper III), the correlation between 
soft X-ray and VLBI radio properties has been discussed. 

In this paper we first give an overview, in § |2l of the sources 
available for the motion analysis. There is a full list of CJF 
source identifications and redshifts, which have not previously 
appeared systematically; we give particulars concerning new 
redshifts in Appendix lAl In §[3]we review the cross-epoch iden- 
tification of components, and discuss the calculation of proper 
motions and their uncertainties. Appendix IB] delves into more 
details about the kinematic model fitting. §|4]discusses the con- 
version to apparent velocities and some further consistency anal- 
ysis. AppendixICjpresents the full table of the kinematic results, 
together with pertinent source and component data. We begin 
the analysis by scrutinizing in §|5] velocity differences between 
components within individual jets, as well as the related topics 



of acceleration and bending. Then, in §|6]we present and analyze 
the velocity statistics of the full sample. § |7] investigates corre- 
lations with other radio parameters, and in particular with the 
observed radio luminosity, while §[8] compares our results with 
other apparent velocity datasets. We give our conclusions in §|9] 

2. The sample 

2. 1. Identifications and redsliifts 

The CJF sample contains 293 radio sources. The properties of 
their optical hosts have not yet been tabulated for all sources in 
previous papers. We present the relevant data here (Table lCTI ) as 
it is required to analyze and interpret the apparent jet component 
motions. 

We have labeled the hosts as either Q(SO), G(alaxy), or 
B(L Lac Object), using one consistent classification from the 
literature, as given by Veron-Cetty & Veron (2003), whenever 
this was available. To classify the remainder of the sources, 
we have inspected optical images, initially the POSS, supple- 
mented later with Palomar 60" and 200" images obtained by 
some of us (Taylor & Vermeulen, unpublished). There are 8 
optical hosts which have never been detected, and are therefore 
U(nclassified); these probably include galaxies at substantial 
redshifts as well as some highly reddened quasars. The method 
of classifying is subject to some uncertainty; in particular, faint 
extended emission, and therefore classification as G, might have 
been missed in some cases. However, these "interlopers" into 
the Q class are expected to be, at most, a minor contaminant of 
that large class. 

We have taken redshifts in preference from our own work 
(Vermeulen & Taylor 1995; Vermeulen et al. 1996; Henstock et 
al. 1995), and we have otherwise used NED to check for the ex- 
istence of a known redshift measurement, which we have traced 
back to the original reference whenever possible for verification 
and maximal accuracy. For 19 sources, our listed redshift (or, 
in one case, the absence of a redshift) constitutes an update to 
previously available information (see Table IC.ll column 3 for 
values). Details for these 19 objects are given in AppendixlAl 

2.2. Sources omitted from ttie motion analysis 

Not all of the CJF sources could be included in the kinematic 
analysis. Table [T] gives a detailed overview of the number 
of sources and jet components that have been used for the 
kinematic analysis, and those that had to be removed from the 
analysis for various reasons. The first column lists the various 
reasons for source elimination, and source name(s) that fell 
victim to each criterion. For each separate reason, column 2 
lists the number of sources affected, column 3 provides the 
cumulative number of rejected sources, and column 4 lists the 
total number of sources that remain available for subsequent 
kinematic analysis. 

The first part of Table [T] summarizes the loss of sources 
during the component model-fit stage, described in Paper I. We 
have found usable information for 266 of the 293 CJF sources 
on the proper motions of 779 jet components. Of the 9 sources 
which could not be used because they are unresolved or have 
a single component in all epochs, as many as 6 are quasars 
at z > 1. Such redshifts are rather common for quasars in the 
sample. The 13 sources in which the components cannot be 
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uniquely identified across epochs do not stand out from the 
sample as a whole in their redshifts or identifications. 

The second part of Table [T] summarizes those sources 
lost from the apparent-velocity analysis (eliminating sources 
without reliable redshifts), and breaks down the number of 
components in all sources remaining in both the proper-motion 
and apparent-velocity analysis. The "quality" of a component is 
discussed in §[3] In Fig. [T] we show the distribution in redshift 
for the 237 sources which finally remain. This sample comprises 
177 quasars, 19 BL Lac objects, and 41 galaxies. 



3. Component identification and proper motion 
fitting 

For all epochs of all sources we first obtained Clean images. We 
then fitted circular Gaussian model components directly to the 
observed complex visibilities using the Levenberg-Marquardt 
non-linear least squares minimization technique (program 
"modelfit" within Difinap, Shepherd 1997). This database is 
presented in Paper I. 

The proper-motion analysis is based on a careful identifica- 
tion of the components across the epochs and multiple checks of 
the resulting motions parameterized in rectangular coordinates 
in the x-y plane. We do realize that "components" are unlikely 
to be sharply defined moving physical entities; rather, they are 
likely to represent those regions in jets which produce a surfeit 
of synchrotron emission as a result of passing shocks and the 
details of the magnetic field configuration. Nevertheless, with 
observing epochs spanning a number of years it is possible in 
most of the sources to trace with confidence the evolution of 
such emitting regions, including their propagation along the jets, 
and this is what, for brevity, is traditionally called component 
motion. When identifying components across epochs in jets 
with multiple features, we followed the assumption that those 
emission regions having the smallest deviations in flux-density, 
core separation, position angle and size between adjacent 
epochs most likely can be identified with the same component. 




Fig. 1. The redshift distribution of the sources included in the 
apparent- velocity analysis (see Table[Tli. The top panel shows all 
237 sources; the lower three panels show separate distributions 
by source class. 
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Table 1. A brief summary of the numbers of sources and components taken along in the subsequent proper-motion and apparent- 
velocity analysis, following through the various stages of elimination. 



(1) 


(2) (3) (4) 


From the complete CJF to the sources with monitored components. 


This Cum. Cum. 
categ. elim. left 


Paper I 




CJF sources: 

Paper 1: Wrong coordinates used for observations: 
0344+405 

Paper 1: Too faint or too extended to image: 

0256+424, 0424+670, 0945+664, 1545+497 

Paper 1: Unresolved or single component in all epochs: 

0621+446 (no z, Q), 0636+680 (z=3.18, Q), 1254+571 (z=0.04217, G) 

1308+471 (no z, VisS by NED), 1342+663 (z=1.351, Q) 

1417+385 (z=1.832, Q), 1638+398 (z=1.66. Q), 1851+488 (z=1.25, Q), 2005+642 (z=1.574, Q) 
Paper 1: components cannot be identified across epochs: 

0205+722 (z=0.895, G): Component identification not sufficiently unambiguous 

0218+357 (z=0.68466, G, BL): Gravitationally lensed 
0402+379 (z=0.055, G): Core identification not unambiguous 
0615+820 (z=0.71, Q): Components not securely identifiable 
0650+371 (z=1.982, Q): Drastic flux-density variability 
0650+453 (z=0.933, Q)i Large position angle changes a.cross epochs 

0954+556 (z=0.895504, Q): Not sufficient epochs due to faintness of source 
1144+402 (z=1.088, Q): No reUable component identification found 
1206+415 (no z, Q) Components not securely identifiable 
1531+722 (z=0.899, G) Blend between core and first jet component 
1800+440 (z=0.663, Q) Resolution of epochs too different 
1839+389 (z=3.095, Q) Jet component too faint 


293 293 
1 1 292 

4 5 288 

9 14 279 

13 27 266 


This Paper 




Proper motion: 

Sources with proper motion measurement(s): 266 

Number of jet components (moving or stationary), all quality classes: 779 
Number of jet components, Ql: 305 (Q: 187, G: 81, B: 26, U: 11) 
Number of jet components, Q2 & Q3: 474 (Q: 299, G: 101, B: 63, U: 11) 
Further kinematic analysis: 
Sources without redshifts: 29 

0018+799 0109+480 0346+800 0633+596 0700+470 0709+619 0716+714 0740+768 

0749+540, 0843+575, 0925+504, 0927+352, 1107+607, 1205+544, 1221+809, 1246+586, 

1250+532, 1312+533, 1322+835, 1333+589, 1357+769, 1656+482, 1746+470, 1747+433, 

1828+399, 1926+611, 2010+723, 2023+760, 2138+389 

Sources with redshifts: 237 (Q: 177, G: 41, B: 19) 

Number of jet components (moving or stationary), all quality classes: 699 

Number of jet components, Ql: 272 (Q: 186, G: 70, B: 16) 

Number of jet components, Q2 & Q3: 427 (Q: 295, G: 93, B: 39) 


29 56 237 
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Cmas) H[gl,t Aacenaidn (mos) Wflht AwMnsTon (moB) 

Contours %: 0.5 0.75 1.13 1.69 2.53 3.8 5.7 8.54 Contours %: 0.8 1.2 1.8 2.7 4.05 6.08 9.11 13.7 Contours %: 2.5 3.5 4.9 6.86 9.6 13.4 18.8 26.4 

Contours %:12.8 19.2 28.8 43.2 64.9 Contours %: 20.5 30.8 46.1 69.2 Contours %: 36.9 51.7 72.3 

Beam F=WHM : 2.29 x 1 .41 (mas) at 1 7.6 ° Beam FWHM : 2.49 x 1 .5 (mas) at -29.4° Beam FWHM : 1 .92 X 1 .44 (mas) at -20.8° 



(a) 



28.08.1995 




1106+380 



08.02.1998 



# 




Right Ascension (mas) 
Contours %: 1 1 .5 2.25 3.38 5.06 7.59 1 1 .4 1 7.1 
Contoufs %: 25.6 3p.4 57.7 86.5 
Beam FWHM : 2.25 x 1 .49 (mas) at -2.64° 



Right AsceriBion (mas) 
Contours %: 1 1 2.25 3.38 5.06 /59 1 1.4 17.1 
Contours %: 25.6 38.4 57.7 86.5 
Bearil FWHM : 2:31 x 1 .37 (mas) at -1 5.4° 

(b) 



21.11.1999 




Right Ascension (mos) 
Contours %: 2 3 4.5 6.75 10.1 15.2 22.8 34.2 51.3 
Contoure %: 76.9 i ' 
Beam FWHM : 2.58 x 1 .71 (mas) at -32.8° 




Right Aacension (mca) ^° ° f ~\° Right Ascension (mas) 

contours %: -3 3 4.5 6.75 40.1 15.2-22.8 34.?51. 3 Contours %: 3 4.5 6.75 10.1 15.2 22.8 34.2 5Y.3 

Contours %:4fr.9 - ! Contours %: -5 5 7.5 1 1.3 16.9 25.3 38 57 85.4 Contours %: 78.9 

Bearri FWHM : 4.81 x 2.5 (rnas) at 1 .07 ° Beann FWHM i 4.81 x 2.5 (riias) at 1 .07 ° 3^,3^ FWHM : 2.98 x 2.72 (masl at -2.03° 



Fig. 2. Images are shown, built up by restoring the model-fitted components, convolved with the clean beam, into the residual image, 
which was made by Fourier transforming the visibility data after first subtracting the model-fitted components in the uv-plane. Over- 
plotted we show symbols to represent the model components, (a) shows the quasar 0711H-356 (all jet components of Ql), (b) the 
galaxy 1 106H-380 (all jet components of Ql), and (c) the quasar 2255-1-416 (CI: Ql; C2,3,4: Q2; C5,6,7; Q3). The components C4, 
C5, and C6 in 2255-H416 split in the last epoch. The identified jet components that can be traced across the epochs are coimected by 
dashed lines. 
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It turned out, however, that not all jet components yield 
proper motion values of equal significance. The uncertainties 
of, and correlations among, the circular Gaussian parameters 
determined in the jet-component model-fitting were propagated 
through the proper-motion analysis (Table IC.lb . but there are 
probably other sources of (systematic) errors as well. A quality 
classification 1, 2, or 3 (high to low) was assigned to every com- 
ponent, based on a careful assessment of the component model 
fits, the images, and the motion in the x-y plane. Single, bright, 
jet components that are clearly separate in all epochs were given 
a quality class 1 . More diffuse, extended, or fainter components, 
and to those in close proximity to others (to the resolution of 
the images) were given a quality class 2. Any component not 
uniquely identifiable at all epochs received quality class 3, which 
also includes a number of components that appeared to merge or 
split, rather than having a single brightness peak at all epochs. 
Such phenomena are to be expected if components represent a 
complex underlying shock geometry, and in the more extreme 
cases have prevented us from deriving a useful motion at all. 
We will denote the quality classification of a component with 
the syntax Q followed by a numeral 1-3 (a "Q" without a trail- 
ing numeral can be used to refer to quasars as a class of sources). 

In Fig. 12] we show examples for jet component identifica- 
tion, quality-class assignment, and motion in three sources: 
0711H-356 (z=1.620, Q), 1106H-380 (z=2.29, G), and 2255H-416 
(z=2.150, Q). The jet components that can be followed across 
the epochs are connected by a dashed line. The component 
labels comprise the letter "C" plus a number increasing with 
ordinal separation from the core (although not illustrated here, 
the prefix "CC" denotes a counter-jet component). The quality 
classification of all components in these sources are Ql, except 
for some in 2255H-416: components C2-C4 are Q2 and C5-C7 
areQ3. 

The position of the reference point (r) is marked by a solid 
black line. We usually call the most compact feature the refer- 
ence point. In a few cases, however, we have chosen a feature 
at one extreme end of the source, even if at some epochs it ap- 
peared somewhat extended, because that could then be ascribed 
to a flare in progress leading to the formation ("ejection") of a 
new feature in the jet. In some obviously two-sided sources we 
chose the most plausible central, core, component to be the ref- 
erence. In others, however, the core was too faint, or unclear, and 
we then used a feature that appears to be the hot spot at the end 
of one of the jets. It is generally assumed that the "core" feature 
at the base of radio jets is stationary, but it has only rarely been 
possibly to investigate this directly with phase referencing (e.g., 
Bartel et al. 1986), and has in fact been called into question in 
some cases (e.g., Guirado et al. 1998). The analysis in this paper 
assumes that the core is stationary to within the errors (but see 
also §|6?T)- 



3. 1. Proper motion calculations 

We calculated the proper motions based on the component 
identifications and model-fit parameters presented in Paper I. 
However, effects based on the dynamic range of the observations 
and source-intrinsic peculiarities complicated this calculation. In 
particular, jet components could be seen to merge or split from 
one epoch to another This can result from comparing observa- 
tions with different dynamic range or angular resolution aris- 
ing from the use of different telescope arrays at different epochs 
(see Table 1 in Paper I). Furthermore, there were some sources 



(a) 



(c) 




Ai 


n=305 




■Wf ■ ■ 




■ iii 


n=474 

Sit'.:. ■ ■ 


■ 





[1 [mas/year] 
(b) 



\i' [mas/year] 
(d) 



Fig. 3. The yu' and fi^' components of the proper motions are 
shown (a) for Ql only, and (b) for Q2 & Q3. The radial and or- 
thogonal proper motion pairs are also shown (c) for Ql 
only, and (d) for Q2 & Q3. "n" is the number of values (compo- 
nents) displayed in this figure. 



in which all jet components could not be traced through all the 
epochs. We investigated several different ways of handling these 
sorts of complications in the estimation of the kinematic param- 
eters for individual jet components. In Appendix iB] we provide 
some details of these procedures but it is appropriate to note a 
few generic points here. 

The kinematic model, for each of the numbered C/CC com- 
ponents, comprises the position {X, Y) at the reference epoch rel- 
ative to the reference point, the proper motion (/i\/i' ), the as- 
sociated uncertainties (ct),^, cr^j) and the correlation matrix. The 
reference epoch is the mean of the observing epochs for that 
source - even if particular components weren't detected in some 
epochs. These kinematic models make no allowance for accel- 
erations or other non-linear motions. The uncertainties in the 
kinematic-model parameters are scaled to a reduced chi-square 
of unity. This procedure helps to allow comparisons among mod- 
els computed using alternate (relative) weighting schemes for 
the uncertainties in the component positions (c/ Appendix [B). 

We also convert the motions for each component from 
an ix,y) to a {radial, orthogonal) representation - both are 
Cartesian frames, but the latter is shifted to each component's 
(X, Y) and rotated such that r points away from the component's 
reference point (i.e., the units of jj" remain mas/yr). The compu- 
tation of the uncertainties in ii/,^") follow from standard error- 
propagation from the original parameters, including correlation 
information. The total proper motion (/i'°') is the quadrature sum 
of the two proper-motion components (in either representation). 

We present the table of kinematic model parameters in 
AppendixICl along with corresponding plots. The complete table 
of model parameters and plots are available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79. 128.5) 
or via htt p://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-HA/, 
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n (qualities 2 & 3)=474 
n (quality 1)=305 



X 0-2 

|4, [mas/year] 



0.4 



(a) 



0.57 
0.49 
0.41 
0.33 
0.25 
0.16 
0.08 





n (qualities 2 & 3)=474 
n (quality 1)=305 



|J, [mas/year] 



0.4 



(b) 



0.57 r 

0.49- 
0.41 - 
0.33- 
0.25- 
0.16- 
0.08- 
0- 



-0.2 



1 



n (qualities 2 & 3)=474 
n (quality 1)=305 



^ 0.2 0.4 
^ [mas/year] 




n (qualities 2 & 3)=474 
n (quality 1 )=305 



0.2 



^ [mas/year] 



0.4 



(C) (d) 

Fig. 4. Histograms showing the distribution of apparent proper motions in orthogonal pairs of values: (a) jt/ ' with (b) and (c) fi'' 
with (d) fi". In each case, Ql components are shown with the solid thick black line, and Q2 & Q3 are in grey scale. The fractional 
distribution was computed separately for each subset, therefore the area under each curve sums to unity. All further histograms in 
this paper have been prepared in this way. 



3.2. Proper motion consistency analysis 

We now show various diagrams which were designed to check 
the consistency of our kinematic modeling and to detect any 
calculation-induced biases. These tests mainly serve as diagnos- 
tics of possible measurement errors but in addition already pro- 
vide some insight into the astrophysics of these sources. 

In Fig. [3] we display the Cartesian motion pairs of the ap- 
parent proper motions: (//*,//') and {j/,^"). The Ql components 
are plotted separately from the Q2 & Q3 components. Fig. |4] 
shows histograms of the same data, and allows a better impres- 
sion of the inner, crowded part of the distributions. As expected, 
there is no preferred direction on the sky. The x,y motions 
form a circular distribution around (0,0), with median apparent 
proper motion values in the Ql components of //*=0.002±0.056 
mas/yr and yU- =0.003 ±0.056 mas/yr The corresponding values 
for the Q2 & Q3 components are yu'=-0.006±0.083 mas/yr and 
//-*'=0.003 ±0.087 mas/yr 

We can also look at the population of components per quad- 
rant in the yu plane to see whether their distribution is consis- 
tent with placement in four "unbiased" bins. The expected stan- 
dard deviation from a binomial distribution with a "probability 



of success" p - IjA can be expressed as crt - -\/np (1 - p) - 

V3«/4. Considering the placement of components per quadrant 
regardless of their proper-motion uncertainties we obtain counts 
of (68,76,81,71), with a mean of 74 components and cr^ - 7.4. 
The Q2 & Q3 components show somewhat more variation in the 
number of components per quadrant: (107,107,123,130), mean 
= 116.75, and cri, - 9.4. The number of components in these 
counts is somewhat lower than the total population in Table [T] 
because components having either fi'' - (7 components) or 
fi>' = (9 components) to the precision of Table IC.ll were ex- 
cluded. If we consider only components falling within a quad- 
rant by at least 1 cr in both jj.^ and yu'', then the distribution of Ql 
components is (29,36,35,31), with a mean of 32.75 components 
and cTfe of 5.0. For the Q2 & Q3 components, the distribution 
is (22,29,26,30), mean 26.75 components and cTb = 4.5. Fig.|5] 
shows the distribution of cr^v. The median of the uncertainties is 
cr^i = 0.022 mas/yr for the Ql components, and cr^jx = 0.049 
mas/yr for the Q2 & Q3 components. These can be compared to 
the standard deviations of the yu^ distribution: 0.078 mas/yr for 
the Ql components and 0.1 15 mas/yr for the Q2 & Q3 compo- 
nents. 
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By contrast, and as expected, the distributions are markedly 
different after conversion to the component specific (radial, or- 
thogonal) frame: clearly, there is a predominance of outward ra- 
dial motions in Fig. [3jc) and (d) (which we did not impose in 
any way as an a priori constraint). There were some compo- 
nents showing inward-directed motion, which will be discussed 
in more detail in § 16.11 

We can also examine the prevalence of clockwise or counter- 
clockwise motion of jet components with respect to the core. 
Fig. m shows a roughly symmetrical distribution of the orthogo- 
nal motions in positive and negative values, as one would expect 
if there is no intrinsic handedness. Following a similar binomial- 
sampling exercise as above, we can quantify this symmetry 
(here, p - 1/2 and thus cri, - ■\/n/2). Considering the placement 
of components in clockwise or counter-clockwise orthogonal 
motion without regard to 07,0 results in counts of (140,158) with 
cr^ = 8.6 for the Ql components and (221,241) with o-j, = 10.7 
for the Q2 & Q3 components. Thus each sub-group has about a 
1-cr offset from the expected symmetry. However, if we consider 
only the components that have values of p° at least cr^„ different 
from 0, then the symmetry is improved: (62,67) with cri, - 5.7 
for the Ql components, and (106,102) with a-j, -1 .2 for the Q2 
& Q3 components. This also shows that over half of the p" for 
all quality classifications aren't significantly different from 0. 

However, the Ql components do have a tighter p" distribu- 
tion than do the Q2 & Q3 components: standard deviations of 
0.04 mas/yr for Ql and 0.10 mas/yr for Q2 & Q3. Along with 
Fig. |5] this supports the appropriateness of the assignment into 
quality classes, and strongly suggests that we concentrate our 
statistical analysis on the Ql components. Nonetheless, the ev- 
idence is that the lower quality components are sampling the 
same motion statistics. We will revisit this issue in SHTI 



4. The database of apparent transverse velocities 

CSapp) 

Apparent transverse velocities (ySapp) have been calculated us- 
ing the 1-year WMAP data (Spergel et al., 2003); the values 
adopted are h = 0.71+° °^; Q.„,h^ = 0.135!° °JJ^; Qtot = 1.02 ±0.2. 
Differences in the apparent velocities due to differences between 
the 1-year and 3-year WMAP parameters are less than 0.1c for 
p - Q.\ masyr"' out to z = 4; this is small with respect to the 
formal uncertainties in the apparent velocities. 

The values calculated for the apparent velocities are listed in 
Table IC.ll As noted in Table [1] components in sources without 
reliable spectroscopic redshifts are omitted from further analy- 
sis. This reduces the number of sources under consideration from 
266 to 237, and the total number of jet components from 779 to 
699. There are 272 Ql components in 150 sources, of which 186 
are in 109 quasars, 16 are in 12 BL Lacs, and 70 are in 29 galax- 
ies. 

Fig. |6] shows the radial and orthogonal apparent velocity 
pairs (j6app,y6°pp). Histograms allow the inner part of the distri- 
butions to be viewed more clearly, and are shown in Fig. |2l 
These diagrams for the apparent velocities are directly equiva- 
lent to Fig. [3] and Fig. |4]for the apparent proper motions. Again, 
we see the predominance of outward radial velocities, but now 
note some inward velocities, and generally much smaller, but 
not always insignificant orthogonal velocity components, with 
equal prevalence of clockwise and anti-clockwise directions. 
Fig. I2] also shows the distribution of the total apparent veloci- 
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the Ql components, panel (b) the Q2 & Q3 components. 



quadratic sum, the total apparent velocity is always positive and 
has a Ricean noise bias against very low values. 

4. 1. Apparent velocity consistency analysis 

Fig.[8]gives an impression of the significance level of the mea- 
sured velocities. For better comparison we removed the "spikes" 



S. Britzen et al.: A Multi-epoch VLBI Survey of the Kinematics of Caltech-Jodrell Bank Flat-Spectrum Sources 



9 



0.25 



-10 



0.16 



n (quality 1)=272 
n (quality 2 & 3)=424 




r 10 
^ [c] 
app 

(a) 



20 



n (quality 1)=272 
n (quality 2 & 3)=425 




„ 10 

^app 
(b) 



20 



n (quality 1 )=272 
n (quality 2 & 3)=424 




* * 10 
plot [c] 

app 
(c) 



rVi n nnni 



and (c) for all 



Fig. 7. The distributions of (a) yS^pp, (b) [5%^, ^,pp 

components. Some of the outliers visible in Fig. |6] fall outside 
the plotted limits. Quality 1 components are shown as a thick 
solid black line, while the distribution of Q2 & Q3 components 
is shown in grey. 



at which result from components having no defined error. The 
median formal velocity uncertainty is crp_^^ - 1.1c for Ql com- 
ponents, and twice as large, crg^^p = 2.0c for Q2 & Q3 com- 
ponents. The apparent velocities in excess of a few c have there- 
fore been measured with high relative significance: most of these 
components really are significantly superluminal. On the other 
hand, a substantial number of lower velocity values have only 
a modest relative significance and most of these motions are ei- 



ther subluminal or, at most, mildly superluminal. Some of the 
low velocities have, however, been measured with fairly high 
relative significance, as can be seen from Fig.|8] Obtaining more 
precise velocities for a larger number of slow components would 
require observations over a much longer time span and/or with 
much higher angular resolution. 

Figs. |6] and [8] show that the distribution of apparent veloc- 
ities for Ql components is tighter than that of the Q2 & Q3 
components, as it was for the proper motions. But now, the 
difference seems to be "just" due to outliers. The orthogonal 
apparent velocity distributions are centered on zero both in the 
Ql components (median yS°pp - -0.1 c), and in the Q2 & Q3 
components (median /3°pp = -0.0 c), as expected. The median 
apparent radial velocity for the Ql components, jSljpp = 2.4 c, 
is actually slightly larger than the median for the Q2 & Q3 
components, y6!^pp = 2.0 c. But this difference is not significant, 
in view of the formal uncertainties in the velocities, discussed 
above. We have decided to carry out all further analysis in this 
paper using only the Ql components, where the measurement 
uncertainties are not only a factor two smaller, but also more 
uniformly defined. 



5. Motion variations witliin sources 

In many sources we have been able to track multiple compo- 
nents and we find that motions within a source are often not 
consistent with a single uniform velocity applicable to all com- 
ponents along a jet. To complement the analysis in section 4 we 
therefore investigated the scatter of apparent velocity, among the 
components composing individual jets, via the median of the 
distribution of standard deviations in yScomp taken over all jets. 
For jets comprising at least three Ql components, these medi- 
ans are 1.97c, 2.02c, and 0.98c for/?'"', /?'', and respectively. 
For jets having at least three components of any quality class, 
these medians become 2.42c, 2.90c, and 1.61c. These scatters 
can be compared to the full p distributions seen in Fig. |7] The 
scatter of apparent velocity in individual jets (comprised of Ql 
components only) is thus twice as great as the formal velocity 
uncertainty of 1 . 1 c and thus the variations are real. 

We now proceed to investigate various intra-jet motion dif- 
ferences. We begin with a study of accelerations along the jet 
(§ 15. Il l and the related question of an apparent dearth of high- 
velocity components at small projected distances from the core 
(§E2]i. We then turn to bending along the jet (§ |5.3l l. We con- 
clude in i^ |5.4| bv discussing ways to obtain a single "representa- 
tive" velocity per source, as required for properly weighted pop- 
ulation studies. 

5. /. Accelerations along individual jets 

Evidence for increasing velocity in pc-scale jets has been found 
before: e.g., M87 (Biretta et al. 1995); 3C84 (Dhawan et al., 
1998) and Cygnus A (Kxichbaum et al., 1998). It is possible that 
the apparent acceleration is not intrinsic but rather results from 
a change in the jet direction or a change in the jet pattern veloc- 
ity unrelated to the jet bulk velocity. In NGC 315, Cotton et al. 
(1999) found an increasing jet velocity in the inner 5 pc from 
the core based on proper motion measurements as well as on the 
sidedness ratio. In FR I radio sources there is observational evi- 
dence that jets decelerate on larger scales than is relevant here, in 
the sense that they are relativistic near the core and become non- 
relativistic within a few kpc (e.g. 3C 449, FeiTeti et al., 1999). 
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Fig. 8. Overview of the formal significance level of the apparent velocities: the ratio of each velocity j6app over its formal uncertainty 
CTg^pp. On the left: scatter diagrams of the formal significance level (ratio of velocity over uncertainty) for all individual Ql compo- 
nents, as a function of their velocity. A substantial number of components with velocities significant at more than ±5cr has been cut 
off. On the right: histograms of the distribution of the formal significances for all Ql components (black line), compared to aU Q2 
& Q3 components (in grey). Significance levels in excess of +5cr are now included. Panels a,b are for the radial velocity, c,d for the 
orthogonal velocity, and e,f for the quadratically summed total velocity. 



Accelerations and decelerations along the lengths of the jet 
as well as transverse velocity profiles are expected from theo- 
retical considerations (e.g., Blandford 2005) and jet simulations 
(e.g., Tiirler et al. 2000). The jet parts that we see are most likely 
shock waves and this introduces additional kinematic complex- 



ity, since the direction and the speed of the emitting plasma be- 
hind a shock front must differ from the kinematic speed of the 
shock. 

With, in general, only 3 epochs per source in our observa- 
tions, we cannot properly investigate acceleration of individual 
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components (d/3/dt). Instead we characterize the acceleration 
along a jet (or counter-jet) in its entirety (d/3/dr) by performing 
a weighted least-squares fit for a linear model to y6app('")- The lin- 
ear coordinate corresponds to the angular r along a jet/counter- 
jet axis via the angular-diameter distance Da - Dm /{I + z). 
As such, this is a global property of the jet, rather than a local 
property of each component. For computing these accelerations, 
we have used {e.g., 1 106H-380 in the plots in Appendix ICt. 
This is a very simple model, and some of the jets with many 
components are clearly not well described by a single accelera- 
tion (e.g., 0633H-734, 1755H-578). In order to be included in this 
analysis, a jet must pass three hurdles: (1) the host source must 
have a redshift (in order to compute the acceleration in units of 
[c/lyr]); (2) it must have at least three components of the speci- 
fied quality class (to ensure the least-squares fit for the accelera- 
tion has at least one degree of freedom); (3) these three compo- 
nents must have been detected in at least three epochs (to ensure 
that the kinematic fit that produced their proper motions [§[TT| 
also had at least one degree of freedom). The source's core was 
not included as a component. Fits were carried out using com- 
ponents of all quality-classes (number of jets per source class: 
quasars=70, galaxies=25, BL Lacs=10; total=105) and within 
Ql components only (Q=14, G=8, B=l; total=23). Galaxies 
seem relatively over-represented in the Ql sample compared to 
quasars, simply because a smaller fraction of quasars have three 
Ql components in a given jet. In the Ql sample as a whole, all 
but two jets have either 3 or 4 components. 

Fig.|9]plots histograms of both the acceleration and velocity- 
intercept (jS'' at r = 0) obtained from the linear fit to y6!|pp(r). 
There is a slight trend towards a positive median acceleration 
and a clear signal of positive velocity-intercept, both statements 
being more significant for the Ql-only sample. The positive 
velocity-intercept condition merely states that the fitted line 
jS!|pp(r) extrapolated back to the core still has a positive value. 
The implication is that, statistically, there is greater acceleration 
in the inner-most regions which the resolution of our observa- 
tions is not able to sample properly. Fig. [TO] shows histograms 
of the acceleration and velocity-intercept for quasars and galax- 
ies from the Ql-only sample. As a class, the quasars show more 
evidence for accelerations having a more positive dfi' /dr and 

5.2. A dearth of high-velocity components close to cores? 

To investigate the tentative evidence for a prevalence of some 
acceleration along resolvable parts of jets in a different way in 
Fig.fTTk we plot the velocities of all Ql components as a func- 
tion of their projected radial distance from the core. There is a 
lack of fast components at distances below a few pc. In some 
sources high velocities have been measured close to the core at a 
wavelength of 7mm. We discuss this in some detail in § 18.21 We 
believe, that the observed dearth here is largely due to at least 
two selection effects. 

The first effect is related to our measurement method. The 
typical time spanned by our series of observations is 3-5 years. 
Any component with a significantly superluminal apparent ve- 
locity will have moved more than a few pc during the moni- 
toring interval, and its average radial position will therefore not 
be within a few pc the core. Only slower components will have 
stayed around near the core long enough for us to have regis- 
tered them there. It is important to stress, that this effect only 
acts to move the faster components to larger radii in our dataset; 
it does not bias the statistics of the velocity values themselves; 
their distribution, (§|6]l, is completely unaffected. 
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panel (b) the velocity intercept at r = 0. The Ql-only jets are 
represented by the thick line and those without quality-class con- 
straint by the shaded area. 



The second effect is more physical in nature. In §|7]we show 
that the apparent velocities are positively correlated with ap- 
parent luminosity and with redshift (two quantities being them- 
selves strongly correlated in flux-density limited samples). The 
smaller radial distances are resolvable only at lower redshifts, 
and therefore in the less luminous sources (usually galaxies, see 
§ |7]l. In these sources we find that larger velocities do not oc- 
cur at all, at any radius from the core. Fig. [12] illustrates that 
the closest components can only be detected in sources with the 
smaller redshifts. A fairly sharp curve delimits the lower-bound 
of r{z) in our sample of components. The dominance of lumi- 
nosity over radial distance in determining the distribution of ve- 
locities is demonstrated by taking successively higher selections 
in luminosity for Fig.fTTb.c.d. 

5.3. Bending along individual jets 

The highest resolution VLBI observations show that in general 
curvature seems to play a more prominent role in the inner- 
most regions of the AGN (e.g. Krichbaum et al. 1994; Britzen 
et al. 2000). More specifically, evidence for precession is found 
in an increasing number of sources (e.g., Britzen et al. 1999a, 
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2005 a&b; Walker et al. 2001; Lister 2001). Jet components in 
AGN exhibit either ballistic motions away from the core or 
curved paths suggestive of streaming motions along a funnel 
(e.g. Homan et al. 2001; Lister 2001). Some of these bent jets 
resemble helical structures in projection, presumably originat- 
ing in precession of the jet nozzle (e.g., Britzen et al. 2005b). 
The precession can have different physical causes for example 
the precession of the accretion disk in binary black holes (e.g., 
Britzen et al. 200 1 a) or a misalignment between the rotation axes 
of the accretion disk and of a single Kerr black hole (Caproni et 
al. 2004). Components can be ejected at different position angles 
(e.g., Britzen et al. 1999a; Qian et al. 2001; Bach et al. 2005), 
initially with ballistic trajectories, and then later on following in- 
trinsically curved paths (Deim, Mutel & Marscher 2000; Stirling 
et al. 2003). 

As opposed to the "global" nature of the acceleration along 
the jet, as discussed above, we can treat bending in a manner 
"local" to each component. Using only components that meet 
the quality criteria, we compute the bending at component C, 
as the angle between the line segments joining components 
C(,_i)-C, and C,— C(,+i). Here, the core is considered, so that 
a bending for the innermost component can be computed. 
However, no bending is associated with the outermost one. The 
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equation for the bending angle in this scheme is: 



fend = cos 



XaXb + yajb 



(1) 
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where 



(i+i) 



(2) 



(similarly for ya and yt,), and an expression for the associ- 
ated uncertainty comes from standard error-propagation (with 
no correlations in the components' X and Y parameters, cfij 13. lb : 
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where the uncertainties for Xa, xt, ya, and yb come from quadra- 
ture sums of the uncertainties in the two components whose dif- 
ference forms these various intermediate variables. The uncer- 
tainty of the core position is taken as 0. The calculated angle 
is always in the range 0-180 degrees, regardless of whether the 
bending is clockwise or counterclockwise for outward motions 
as exemplified in 2255H-416 (Fig.|2]i by the component sequences 
(C2-C3-C4) and (C4-C5-C6), respectively. 

This bending is characterized above as a property of the in- 
dividual components, but it of course also depends on each com- 
ponent's neighbors that meet the quality criterion. The Ql-only 
sub-set can be expected to have fewer components, spaced far- 
ther apart along the jet thus producing smaller bending angles as 
defined above. This seems borne out in the two panels of Fig.[T3l 
especially for the quasars. 

The size and quality of our data set allows us to take the 
bending analysis one step further Fig.[T4]shows the relation be- 
tween lavT^ {J}" I P' ) and the bending angle 0bnd for Ql-only com- 
ponents and for all components. The former quantity relates to 
the angle a component's velocity makes with its local radial di- 
rection: 0° implies purely outward radial motion, +180° purely 
inward radial motion, and +90° purely circumferential motion 
around the core. Apparent circumferential motion is rare in that 
there are very few Ql components (~4% of the components 
from galaxies and <2% of those from quasars) within 15° of 
+90°. When all quality classes are considered, more components 
show apparent circumferential motion, 5% from quasars, ~7% 
from galaxies, and preferentially more from BL Lacs (~21%). 



Fig. 13. (a) Distributions of measured local angles of bending 
(see text for details). Quasars: thick solid line, galaxies: thin 
solid line, BL Lac Objects: grey-scale. Panel (a) shows results 
when using all components, panel (b) when using only Ql com- 
ponents. 



However, there is no clear correlation between large bending an- 
gles 0bnd and the angle a component's velocity makes with the 
local radial direction. 

5.4. Can a representative velocity for eacli source be 
determined ? 

As discussed in the introduction to this section, there are many 
sources in which we have been able to track multiple compo- 
nents, and in any single source these components often have sig- 
nificantly different apparent velocities. For population studies, 
a single velocity per source must be derived to give each source 
equal weight, for example when comparing the velocity statistics 
of galaxies, where there are often many components per source, 
to quasars, which often have fewer components in their jets. 

The selection from our sample of either the brightest or the 
closest Ql component per source yields a distribution of veloc- 
ities which cannot be statistically distinguished from that of our 
full sample of Ql components. It is important to point out, how- 
ever, that the distribution of brightest Ql component velocities 
terminates at ~ 18 c. Histograms of the relevant distributions are 
shown in Fig. [15] Note that the selected sub-samples of the 150 
brightest or 150 closest Ql components are part of our full sam- 
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Fig. 14. X.?ivr^(ff IP'') as a function of the bending angle 6'bnd for 
the Ql components (a) and all components (b). Quasars are de- 
noted by squares, galaxies by triangles, and BL Lacs by dia- 
monds. 



pie of all 272 Ql components, and the histograms are far from in- 
dependent; the brightest and closest sub-samples overlap for all 
but 23 sources between the brightest and the closest sub-sample. 
Despite the similarities in Fig. [15] we favor the use of the bright- 
est rather than the closest component. With limited resolution the 
position of components close to the (usually time-variable) core 
may be affected by deconvolution problems. We have not given 
Ql labels to components where we suspected "core-blending" 
in our sample. There are also selection effects which increase 
the fraction of slow components at the smallest radial distances 
(<5 pc, see ii l5.2l l. While the full sample of closest Ql compo- 
nents spans a large range of radial distances, mostly above a few 
pc, and so should not be much affected, the use of the brightest 
components might nevertheless be preferable partly because of 
this effect. Of course, if jets bend, or vary in velocity with time 
or along their length, then selecting the brightest components 
would tend to bias towards features which are the most strongly 
Doppler boosted into our line-of-sight. These do not necessarily 
have either the fastest or any other characteristic apparent ve- 
locity (for reference, see the discussion in Vermeulen & Cohen 
1994). 

We have thus shown that taking the velocity of the brightest 
components - this means one velocity per source - as representa- 
tive may be adequate for population studies, knowing that their 
distribution matches that of all Ql components. 



6. The distribution of apparent radiai veiocities in 
thie sampie 

The radial velocities in Fig. [16] (an expanded version of Fig. [6]) 
and Fig. [7] show a fairly broad distribution of velocities, tail- 
ing off towards ~ 15 c, with small sprinkling of radial velocities 
extending up ~ 30 c. The distribution also extends, seemingly 
continuously, to somewhat negative (i.e., apparently inward) ra- 
dial velocities, but these have tapered off at about -5 c. We find 
a pronounced narrow peak of subluminal or, at most, mildly su- 
perluminal sources/components. 

To obtain such a broad distribution of apparent veloci- 
ties, predictions from relativistic jet models (e.g., Vermeulen & 
Cohen 1994; Lister & Marscher 1997; Cohen et al. 2007) re- 
quires a broad distribution of the Lorentz y factors of the moving 
features. 

6. 1. Negative radial velocities 

Of the 272 Ql components in our sample, 56 in 42 differ- 
ent sources (26 of 186 quasars and 16 of 70 galaxies) have a 
(slightly) negative radial velocity. But only 15 of the negative 
Ql component velocities (in 8 quasars and 7 galaxies) have a 
formal significance better than 2cr; the median velocity of those 
is -2.6 c (see also Fig.[8]l. We believe it is likely that some of the 
observed negative radial velocities really represent features mov- 
ing towards the core, for example due to reconfinement shocks 
(e.g., Britzen et al. 2005a). Apparent negative velocities in jets 
have also been discussed by Wehrle et al., (2001), Kellermann 
et al. (2004), and Finer et al. (2007). A theoretical model of ap- 
parently backward-moving knots has been presented by Istomin 
& Pariev (1996) and Gomez et al. (1997). In many cases, how- 
ever, the negative velocities either are simply the result of our 
measurement procedure, or have a geometrical rather than a di- 
rectly physical explanation in the source. We first identify three 
possibilities relating to our adopted measurement procedure: 

1) CSOs (Compact Symmetric Object): there are a number 
of well-known CSOs and similar two-sided sources, often galax- 
ies, where, rather than to the quite inconspicuous or absent core 
at 5 GHz, we have referenced to one of the hot spots, often at 
one extreme end of the structure. Components in the jet moving 
towards that hot spot will then have the sign of their velocity re- 
versed; the hot spot is likely to be almost stationary or at least 
moving much more slowly than components in the jet. Examples 
of such sources include the quasar 2021H-614 and the galaxies 
1946H-708 and 2352H-495. As expected, the absolute values of 
the velocities in the two-sided sources are small. 

2) Incorrect cores: even in some one-sided sources we may 
have, adopted a prominent compact jet component as the refer- 
ence, particularly if it seemed to be (nearly) stationary. Many 
such stationary features exist in jets even amongst moving com- 
ponents, for example as a result of shocks or bends in jets; see 
also 16.21 With the choice of such a reference feature, any com- 
ponent between it and the true core will have the sign of its ve- 
locity reversed, compared to the regular convention where posi- 
tive radial velocity implies motion away from the core. 

3) Blended cores: we also suspect that in some sources there 
was a gradual centroid displacement of the reference feature dur- 
ing the time spanned by our observations, as a result of blending 
between the true core and a component which was newly emerg- 
ing and/or becoming relatively brighter This can also lead to ap- 
parently negative velocities for components which in reality are 
stationary or slowly moving out. 
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Fig. 15. The jS^pp, yS^pp, and j6!,pp distributions obtained with different selections of a single component velocity per source, shown 



app' 



app 



in grey: (a-c) the brightest Ql component in each of 150 sources, and (d-f) the Ql component closest to the core in each of 150 
sources (these sub-samples overlap in all but 23 sources). For comparison, the distribution for all 272 Ql components in the fuU 
sample is shown with a solid black line in all panels. 



Geometric explanations are also possible. Relativistic jets 
are preferentially seen at small angles to the line-of-sight as a re- 
sult of Doppler favoritism. Thus, even slight intrinsic bends can 
create sections of the jet in which components appear to move 
back towards the core in projection. 

6.2. Subluminal and stationary components 

The uncertainties of the motions measured in our sample are 
such that we cannot discriminate well between "stationary" com- 
ponents and components moving at a subluminal or at most 
mildly superluminal velocity. It can be anticipated that samples 
selected on beamed emission (as the CJF was designed to be) 
should show relatively few, if any, subluminal motions (e.g., 
Vermeulen & Cohen 1994). On the contrary, however, we find 
that there is a sharp peak at low velocities. The prevalence of 
subluminal or even stationary components is enhanced even over 
the broad velocity distribution which, by its broadness, is itself 
already difficult to accommodate within simple beaming mod- 
els. In the sample of Ql components, a subluminal radial veloc- 
ity has been measured for 42 of 272 components in 36 sources: 
26 components in 25 quasars, 14 components in 9 galaxies, and 
3 components in 2 BL Lac objects. For Ql & Q2 & Q3 26 sub- 
luminal components in 20 sources have a relative significance of 
2cr or better. 

Slow components occur relatively often, but by no means 
exclusively, in galaxies. However, this cannot be taken as direct 



evidence for standard unification models, in which the relativis- 
tic jets in galaxies are supposed to be viewed at larger angles to 
the line-of-sight than in quasars. For example the CJF sample 
contains a number of well-known CSOs which are always iden- 
tified with galaxies; the prototype is 2352H-495 (e.g., Wilkinson 
et al. 1994; Taylor et al. 1996b; Readhead et al. 1996). It has 
been firmly established that CSOs are not dominated by beamed 
emission, and that components in their jets are typically sublumi- 
nal; speeds of 0.2-0.4 c are common (e.g., Polatidis & Conway 
2003). 

We find that slow jet components often coexist with super- 
luminal components, and can occur in sources at any observed 
luminosity. Such (near-)stationary features probably have a geo- 
metrical explanation in some cases, and a truly physical cause in 
others. Some slow components are likely to represent sites of en- 
hanced Doppler beaming, where a curved jet points most closely 
along the line of sight, as in the source 4C39.25 (Alberdi et al. 
2000). Others are probably associated with a stationary shock 
in the jet, for example a standing recollimation shock caused by 
pressure imbalances at the boundary between the jet fluid and the 
external medium, where the energy density enhancement pro- 
duced downstream can give rise to stationary radio knots (e.g., 
Marti & Miiller 2003). If a jet bends abruptly, whether because 
of internal instabilities, or because the jet is influenced by gradi- 
ents in the external medium, perhaps even to the extent of being 
deflected, it is also likely to develop some stationary features as- 
sociated with shocks; see for example the work by Gomez et al. 
(2001), Gomez (2005) and Agudo et al. (2001). A proper charac- 
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Fig. 16. The radial and orthogonal values (fi'^pp, /3°pp) of the ap- 
parent velocities for the Ql components, on an expanded scale 
with respect to Fig.|6l 



terization of "stationary" components requires a denser sampling 
of the observations in time than is available for CJF, as demon- 
strated by the oscillations in the position of a "stationary" bright 
jet component in 1803H-784, uncovered by Britzen et al. (2005a). 



7. Correlations with iuminosity and otiier source 
parameters 

In this section we first demonstrate that in the flux-density lim- 
ited CJF sample there is a correlation between the distribution 
of apparent velocities and the observed radio source luminosi- 
ties (i) |7.1t . In 17.21 we show the impact of this correlation on 
the overall velocity distributions of CJF quasars, galaxies, and 
BL Lac objects, and briefly discuss this, together with prop- 
erties such as spectral index and core dominance, in the light 
of orientation-unification models. The correlation with observed 
luminosity must also be taken into account when comparing 
samples with different selection criteria (see §|8ll. 

7.1. A correlation with tine observed 5 GHz core 
luminosity 

With the large, homogeneously analyzed flux-density limited 
CJF sample available, the /i-z and y6!,pp-z diagrams, shown in 
Fig. [T7] reveal rather different features than their predecessors 
(e.g., Cohen et al. 1977, 1988), which were based on smaller, 
heterogeneous samples. The left-hand panels of Fig. [Tt] contain 
all Ql components. Each source can have different numbers of 
components, so in order to give all sources (i.e., all independent 
redshifts) equal weight, we also include on the right-hand side 
of Fig. [17] and subsequent figures diagrams which include only 
the brightest Ql component per source. 

The CJF sample shows a well-defined upper envelope in 
yS!j°p-z, particularly sharply visible for the brightest Ql compo- 
nents. However, that upper limit, far from being constant, is a 
strong function of redshift. The CJF sample has a strong corre- 
lation between redshift and source luminosity, inevitable for a 
flux-density limited survey. This is shown in Fig. [18] We have 
derived the apparent radio luminosity (L = F,„4nDj) from the 
observed 5 GHz radio core flux density F,,,, measured through 
model-fitting in the last available observing epoch; procedures 
and results are given in Paper 1. The luminosity distance is com- 
puted as Di = £>m(1 + z), using the Dm from in § |4] Fig. [19] 



shows file correlation of apparent velocities with the observed 5 
GHz luminosities, for all CJF Ql components on the left, and for 
the brightest Ql components on the right. Globally, the same dis- 
tribution is seen in /3fp^-L as for the diagrams with redshift. Th e 
upper envelope is agam particularly well-defined for the bright- 
est Ql components. It runs in the sense that at low luminosities, 
there are no fast motions at all, while the largest velocity encoun- 
tered increases with lumino sity. However, even though clearly 
subject to small number statistics, we note that the highest veloc- 
ities may not occur at the very highest luminosities; a turnover 
of the envelope is part of "aspect curves" as summarized below. 
Despite what seems to be a rather clear envelope, the motion val- 
ues do not crowd up strongly against it. Conversely, they seem to 
be a bit sparser close to the limit, while the rest of the "allowed" 
velocity range at any luminosity seems covered fairly uniformly 
almost down to zero motion. 

The shape of the Pfpp-L has meanwhile been confirmed re- 
markably consistently by Cohen et al. (2007) in the 2cm VLBA 
survey, which has an overlap of only 24 sources with the CJF. 
Cohen et al. (2007) give an extensive discussion of beaming 
models which can reproduce the global shape of the upper enve- 
lope remarkably well with so-called "aspect curves", that trace 
how the observed luminosity and speed would vary as a given 
source is turned from having its jet pointed directly at the ob- 
server (maximal luminosity, no speed), through small angles 
(slightly reduced luminosity, maximal speed), to large angles 
(strongly reduced luminosity, low speed). However, Cohen et al. 
(2007) have concluded that, in reality, to explain the detailed dis- 
tribution of velocities as a function of observed luminosity, many 
different aspect curves are required. 

Since the observed "high frequency" (i.e. > 5 GHz) lumi- 
nosity is a quantity which is physically closely tied up with the 
relativistic jets that also display the moving components, we 
are confident that luminosity, rather than redshift, is the vari- 
able which truly correlates with velocity. A secondary argument 
in favor of this assertion is that much of the gradient in ve- 
locity as a function of redshift takes place over cosmologically 
fairly late times, certainly compared to "the great quasar era" 
(1.5<z<3). Breaking the luminosity -redshift degeneracy directly 
requires a restricted range of luminosities with a sufficient num- 
ber of sources spread over a reasonably large range in redshift, 
or an orthogonal cut restricted in redshift. Fig.|20]shows the best 
such selection afforded by the CJF sample: these are /3fpp-z dia- 
grams in which only sources in the restricted range of luminosi- 
ties 10^^-10^"' ergs ' Hz are included. There is no longer a 
definite sloping upper envelope to the velocities as a function of 
redshift which supports the idea that, indeed, the primary corre- 
lation is with observed luminosity. 
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Fig. 17. Motions as a function of redshift. On the left are all Ql components, on the right the brightest Ql component per source: 
(a,b) the well-known "ju-z" diagram of apparent proper motions (//'°') against redshift; (c,d) the same relationship with a logarithmic 
redshift axis; (e,f) the proper motions converted to speed OS^pp) as a function of redshift (logarithmic axis). The quasars are indicated 
by filled black squares, the BL Lac objects by open circles, and the galaxies by striped triangles. 
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Fig. 18. The apparent luminosity, calculated based on the 5 GHz 
core flux density, as function of the source redshift for all sources 
with a Ql jet component. The quasars are indicated by filled 
black squares, the BL Lac objects by open circles, and the galax- 
ies by striped triangles. 



7.2. Other correlations 

Apparent superluminal velocity studies have often indicated dif- 
ferences in the statistics between objects classified as quasars, 
galaxies, and BL Lacs (e.g., Gabuzda et al. 1994; Vermeulen 
1995; Jorstad et al. 2001 ; Finer et al. 2007). We show histograms 
of the distributions in these different classes for the CJF survey 
in Fig. |2T1 The CJF galaxies clearly have slower apparent jet 
component velocities than the quasars, on average. We have ver- 
ified that this conclusion holds when the known CSOs are left 
out. 

However, we believe this trend may not be intrinsic to the 
classification of the host object. Fig.[T9lshows an increase in ySl,°p 
with observed luminosity for all classes of sources; the same is 
true as a function of redshift (Fig. [TTI l. Over the range where 
the quasars and galaxies overlap in redshift and luminosity, their 
motion distributions cannot be distinguished. Indeed, superlumi- 
nal velocities up to > 10 c exist in some galaxies. Classification 
of the optical host type (we have taken these from the literature 
in many cases) is often an uncertain undertaking: Much further 
work is needed to clarify how the ratio between stellar and non- 
stellar optical continuum, radio continuum, and jet component 
velocities are inter-related. 

There are too few BL Lac objects with a known redshift in 
the CJF sample to be able to study their velocity distribution in 
comparison to quasars or galaxies. In the past, different studies 
have reached opposite conclusions on whether BL Lac objects 
typically have slow or fast jet components (e.g., Gabuzda et al. 
1994; Wehrle et al. 1992; GhiselHni et al. 1993; Vermeulen 1995; 
Jorstad et al. 2001). We suspect this can arise when different 
definitions of BL Lac objects are used. 

We have also investigated whether the velocities in CJF 
sources depend on their 1 .4-5 GHz radio spectral index, or on 
the ratio at 5 GHz between the VLBl-core flux density and the 
single-dish (i.e. total) flux density ("Core-dominance") . In the 
context of quasar-galaxy orientation unification models, these 
parameters are sometimes taken as proxies for the amount of 
beaming, or for the jet angle to the line-of-sight. But for the CJF 
sample, neither of these properties is correlated with the compo- 
nent velocities. 



8. Comparisons between motion surveys 

Compared to the CJF survey, previous compilations of apparent 
velocities contained far fewer sources, and usually were rather 
more heterogeneous, both in the sample composition, and in the 
observing setup. In contrast to our CJF results, these older lists 
typically contained a large fraction of highly superluminal veloc- 
ities, many in excess of 10c, and ranging up to 40c and more (af- 
ter conversion to current cosmology). There may also have been 
a preference to study highly variable objects (e.g., Wehrle et al. 
1992), which may tend to select for high Doppler factors (e.g. 
Homan et al. 2006). The earlier observations typically involved 
sources which are several times brighter than most CJF sources. 
Those sources at significant redshifts are often amongst the most 
luminous ones around. The correlation discussed in 17.1 [ implies 
a significantly enhanced chance to have a fast velocity at the high 
end of the luminosity distribution. 



8. 1. Comparison witli tlie VLB A 2cm Survey and tlie 
RRFID Survey 

Finer et al. (2007) have recently published the kinematics of a 
significant sample of 77 AGN taken from 8 GHz VLBA images 
in the RRFID (Radio Reference Frame Image Database). These 
data are based on the first 5 years of the database (1994-1998). 
The overlap between the CJF and this stage of the RRFID survey 
is 17 sources. A larger sample of velocities has been published 
for 110 sources from the VLBA 2cm survey (Kellermann et al. 
1998, 2004). While this is still only about half as many sources 
as in CJF, the observing and data analysis were carried out ho- 
mogeneously, and a large number of observing epochs were ob- 
tained for each source. The overlap between the samples is 24 
sources with measured motions. 

The VLBA 2cm survey has brighter sources (limit 1 .5 Jy on 
the extrapolated 2cm flux density), no spectral index limit, and 
covers a much wider area of sky compared with the CJF sample. 
A cursory inspection shows that the velocities in these common 
sources often differ significantly between the two surveys. But 
the observations often do not pertain to the same components 
because the measurements often did not span the same time pe- 
riod and because the two surveys differ in spatial resolution by 
about a factor of three. Disentangling the situation in the individ- 
ual sources is beyond the scope of the present work, and requires 
concurrent multi-frequency observations in the future. We can, 
however, compare the velocity statistics of the two full samples. 

We compare the histograms of the apparent velocity distri- 
butions i.e. Fig. |7] in this paper, and Fig. 5 in Kellermann et 
al. (2004). In direct analogy to the distribution of velocities dis- 
cussed in for the CJF distribution, we find that the VLBA 2cm 
survey distribution shows evidence for two populations: one with 
slow or at most mildly relativistic components, the other with a 
broad range of velocities extending to highly superluminal val- 
ues. Indeed, Kellermann et al. (2004) have already proposed that 
the "stationary" components may well represent a different phe- 
nomenon. The measured apparent speed distribution for the 94 
best-measured components (Fig. 4) by Finer et al. (2007) also 
shows a peak at low apparent speeds together with a tail extend- 
ing out to about 30c. 

Returning to the VLBA 2cm survey, the distribution of the 
superluminal velocities appears to be flat out to almost 10c, 
whereas the distribution in CJF already turns down beyond 4-5c. 
Likewise, while velocities in the range 15-25c are very rare in 
CJF, they are more common in the VLBA 2cm survey. Thus, the 
mean velocity in the VLBA 2cm survey is substantially higher 
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Fig. 19. The total apparent velocities as a function of the observed source core luminosity at 5 GHz: (a) for all Ql components; 
(b) for the brightest Ql component in each source. The quasars are indicated by filled black squares, the BL Lac objects by open 
circles, and the galaxies by striped triangles. 
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Fig. 20. The total apparent velocities as a function of redshift, restricted to sources for which the observed core luminosity is in 
the range lO-'^-lO^^ erg s"' Hz ^ (a) shows all Ql components, (b) only the brightest Ql component per source. The quasars are 
indicated by filled black squares, the BL Lac objects by open circles, and the galaxies by striped triangles. 



than in the CJF survey; discounting the probably separate pop- 
ulation of "stationary" components, the ratio is roughly a factor 
of 1.5-2. 

The correlation discussed in ij 17.1 [ between the apparent ve- 
locities and the observed luminosities, is also prominent in the 
VLBA 2cm survey; compare Fig. [T9]to Fig. 6 in Kellermann et 
al. (2004). There is a larger proportion of highly superluminal 
velocities in the VLBA 2 cm survey than in the CJF sample. 
However, the luminosities in the VLBA 2cm survey are slightly 
larger, as one might expect given that the sources are brighter. 

8.2. Increasing apparent velocity with increasing 
observing frequency 

In the context of increasing apparent velocity with observing fre- 
quency a survey including results at 7mm has been published by 
Jorstad et al. (2001). Both the mean and the largest velocities 
are higher than in the VLBA 2cm survey and we have already 
noted that the VLBA 2cm survey apparent velocities seem to be 
faster than the CJF-results. However, while Jorstad et al. (2001) 
have published an impressive body of work, their sample size is 
comparatively modest. And with its focus on sources detected in 
gamma-rays, it is likely to be less representative of a complete 
flux-density limited survey. It does provide evidence that sources 



prone to being detected in gamma-rays often show apparent su- 
perluminal motion (e.g.. Bower et al. 1997; Barthel et al. 1995. 
There are plausible beaming models which can explain the clear 
trend for higher velocities to be seen at higher frequencies. For 
example, if individual jets contain matter moving with a range of 
velocities, perhaps in a fast spine and slower sheath configura- 
tion (e.g., "two-fluid model": Pelletier & Roland 1998; Laing & 
Bridle 2002), then optical depth effects could reveal faster mate- 
rial at higher frequencies. Alternatively the link between velocity 
and resolution could be taken as evidence that higher velocities 
occur closer to the central engine indicating that some decelera- 
tion already takes place on parsec scales. This, however, conflicts 
with the very tentative evidence for a greater prevalence of ac- 
celeration rather than deceleration along the jets seen in our CJF 
data. Overall, then, while the velocity-frequency trend is clear 
much further work is needed to understand the details of the ef- 
fect and how it impinges on the theory of relativistic jets. 

9. Conclusions 

In this paper we have described the statistical analysis of jet com- 
ponent motion in 237 CJF-sources. We have described in detail 
the necessary tests that have been performed to ensure that our 
treatment is free of any calculation-induced bias. We have inves- 
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Fig. 21.j6!,°p for the different source classes based on the bright- 
est sub-sample. The quasars are indicated by a thick solid black 
line, the BL Lac objects in grey, and the galaxies with a thin 
black line. 

tigated several different ways of handling various sorts of com- 
plications in the estimation of the kinematic parameters for indi- 
vidual jet components, for example splitting and merging com- 
ponents mergers. 

The program we developed allows us, from the computation 
of 2-D kinematic models, to estimate both radial motions 
and orthogonal motions. The overwhelming predominance of 
outward radial velocities in AGN is confirmed, as well as the 
dominance of significant apparent superluminal motions. Our 
most important results are briefly summarized as follows: 

a) Observational results from CJF alone 



Bending along individual jets: 

Mild degrees of bending are common but we find very 
few Ql components undergoing apparently large bends i.e. 
within 15° of +90°. Only ~4% of the components from 
galaxies and <2% of those from quasars show such bends. 

Representative velocities: 

Despite the subtleties of the motions within a given jet, a 
single velocity per source is required for population studies. 
We find that the sample of the brightest Ql components (one 
per source) yields a distribution of velocities which cannot 
be statistically distinguished from that of our fuU sample of 
Ql components; this has therefore been used as the velocity 
sample for our various correlation studies. 

Radial velocity distribution: 

The distribution of radial velocities shows a broad distri- 
bution of velocities (apparent velocities up to 30 c) with 
a pronounced narrow peak of slow apparent velocities 
superimposed. These subluminal and stationary components 
occur relatively often: 42 of 272 components of the Ql 
components reveal such "slow" radial velocities which 
may need to be explained in a different manner to the 
superluminal motions. For the latter the broad distribution 
of apparent velocities requires a broad initial distribution of 
the Lorentz factors of the moving features. 

Negative radial velocities: 

Negative or "backwards" superluminal motions are seen 
and in 15 cases these are significant. In some of these cases 
the phenomenon could be ascribed to geometrical effects 
in curved jets seen close to the line of sight or to actual 
backward motions of shocks. However we suspect that in 
the majority of cases the apparent backwards motions arise 
from problems with our choice of reference feature. 



- Motion within sources: 

In general, we find a spectrum of motion values in a given 
source: motions are not consistent with a single uniform 
velocity applicable to all components along a jet. The 
median of the distribution of standard deviation taken for 
jets comprising at least three Ql components are 1.97c 
(yS""), 2.02c (ySO, and 0.98c (/3°). 

- Accelerations along individual jets: 

Although three epochs are generally not sufficient to prop- 
erly investigate acceleration of individual components along 
the jet, we can characterize the acceleration along a jet (or 
counter-jet) in its entirety. For the Ql-only sample we find 
a slight trend towards a positive outward acceleration. We 
also adduce some evidence for greater acceleration in the 
inner-most regions but these observations do not have the 
resolution to test this properly. The quasars reveal a more 
positive dp'' Idr and fi'' (r = 0). 

- A lack of fast components at distances below a few pc: 
This effect might arise simply from the limited time 
sampling of the motions; this could easily be tested in a sub- 
sample observed with smaller time intervals. Alternatively 
the effect could be physical in nature relating to the fact 
that the smaller physical distances are resolvable only at the 
lower redshifts and thus in the less luminous sources, where 
we do not find large velocities. 



b) Correlations with other source parameters 



- Correlation with the observed 5 GHz luminosity: 

We have been able to break the degeneracy between redshift 
and luminosity and find a strong correlation between the 
5 GHz luminosity and apparent velocity. A clear outcome 
of this analysis is that at the high end of the luminosity 
distribution the chance to have a fast velocity is significantly 
enhanced. The correlation with luminosity may contribute to 
the lack of fast components at the smaUest core separations. 

- Apparent velocity with regard to the source class: 

The CJF galaxies clearly show slower apparent jet com- 
ponent velocities than the quasars, on average. However 
this difference may well be due to the velocity-luminosity 
correlation since over the range where quasars and galaxies 
overlap in redshift and luminosity, their motion distributions 
cannot be distinguished. The number of BL Lac objects with 
a known redshift in the CJF sample is too small to aUow a 
meaningful comparison with quasars or galaxies to be made. 

- No correlation with (1.4-5 GHz) radio spectral index or the 
core-dominance parameter: 

Neither of these parameters shows any correlation with the 
component velocities. 



S. Britzen et al.: A Multi-epoch VLBI Survey of the Kinematics of Caltech-Jodrell Bank Flat-Spectrum Sources 



21 



- Dependence of the average apparent velocities on the 
observational frequency of the motion survey: 
A comparison of the results of the CJF with previous 
compilations of apparent velocities is hampered by the fact 
that they contained many fewer sources. The dependence 
of velocity on the luminosity additionally complicates a 
comparison. Nevertheless the mean velocity in the VLBA 
2cm survey is substantially higher than in the CJF survey 
- the ratio could be roughly a factor of 1.5-2. This trend 
of increasing apparent velocity with increasing observing 
frequency is supported by even higher velocities found 
by Jorstad et al. (2001) in a VLBI survey observed at 
frequencies up to 43 GHz. Simultaneous multi-frequency 
studies of well-chosen sub-samples are needed to clarify 
this phenomenon. 
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Appendix A: New Redshift Information 

For 19 sources, our redshift determinations constitute an update 
to previously published values (see Table IC.ll column 3 for val- 
ues). Details for these 19 sources are as follows: 

0018+729 

An unpublished Palomar 200" Hale Telescope observation by 
Vermeulen & Readhead on 1994 August 15 does not confirm the 
presence of the [O Il]/13727 line on the basis of which Snellen et 
al. (1996) reported a redshift z=0.821. 

0340+362 

The probable emission line redshift published in Vermeulen & 
Taylor (1995) was confirmed in an unpublished Palomar 200" 
Hale Telescope observation by Vermeulen & Readhead on 1996 
January 16, which revealed the C Ill]/I1909 emission line in ad- 
dition to some of the emission lines already tentatively identified 
in Vermeulen & Taylor (1995); the final redshift measurement is 
z=1.485±0.002. 

0346+800 

This source still has no secure emission line redshift of 
which we are aware. However, an unpublished Palomar 200" 
Hale Telescope observation by Vermeulen & Readhead on 
1996 January 14 shows a probable narrow emission line near 
7550 and a possible faint broad emission feature near 5500 . 
If interpreted as [O Il]/13727 and Mg II/J2798, respectively, these 
would tentatively suggest a redshift near z=1.03. 



The probable redshift, z=l .95, already listed in Kellermann et al. 
(1998), is based on an observation at the Keck I Telescope by 
Taylor, which shows irregularly shaped emission features con- 
sistent with C IV^1549 and C III] 1909. 

0604+728 

The emission line redshift, z=0. 986+0. 001, is based on an 
unpublished Palomar 200" Hale Telescope observation by 
Vermeulen & Readhead on 1996 January 14. It shows narrow 
[O Il]/13727 and [Ne V]/13426 emission lines, and more tentative 
H5/14102 and [Nelll]^3869 emission Hnes. Mg 11^2798 is also 
seen in emission, but has an associated absorption line compo- 
nent as well; this associated absorption system is also visible in 
Fe IU2586, 2600 and i2374,2382 hnes. 

0738+491 

The probable emission line redshift, z=2.32, is based on an 
unpublished Palomar 200" Hale Telescope observation by 
Vermeulen & Readhead on 1996 January 14. It shows a sharp 
emission line near 4040 , which is probably Lyci'/I1216, and a 
broad, irregular emission feature near 5 140 that probably corre- 
sponds to ClV/il549. 

0800+618 

The emission line redshift, z=3.044±0.002, is based on an 
unpublished Palomar 200" Hale Telescope observation by 
Vermeulen & Readhead on 1996 January 14, which shows 
LyaAUie, Si IV/O IV^1400, and C IV^1549 in emission. There 
is also evidence for emission lines corresponding to Lyy6/11026 
and Clll]/11909. Furthermore, there is separate, deep, possi- 
ble damped, absorption line system at Zabs=2.963, in which 
at least the following lines ai-e visible: ClV/11549, Si 11/11527, 
CIU1334, OU1302, Sill/11216, LyQ'.11216, and Lyy3.11026. 

0942+468 

The emission line redshift, z=0.639±0.002, is secure. It is based 
on naiTow [Oll]/13727 and [O Ill]/14959, 5007 emission lines 
in an unpublished Palomar 200" Hale Telescope observation by 
Vermeulen & Readhead on 1996 January 16. Engels et al. (1998) 
have published another value, z=0.993, but this was based on as- 
suming that a single feature centered at 5577 (where there is a 
very prominent night sky line) would be Mg 11/12798. 

1125+596 

The emission line redshift, z=l. 799+0. 003, is based on unpub- 
lished Palomar 200" Hale Telescope observations by Vermeulen 
& Readhead on 1996 January 14 and 16, which show an 
unambiguous combination of emission lines: Si IV/O lV/11400, 
C lV/11549, and C Ill]/I1909. There is also evidence for an emis- 
sion line corresponding to Mg 11/12798. The resonance lines have 
an irregular profile due to several unresolved associated absorp- 
tion systems which cut into the broad emission lines. 

1300+580 

Unpublished observations by Vermeulen, Taylor, & Readhead 
at the Palomar 200" Hale Telescope on 1994 May 15 and 1996 
January 16, and at the Keck I Telescope on 1996 June 18, to- 
gether show the emission line redshift to be z=l. 088+0. 003, 
based on Mg 11/12798 and [O Il]/13727 emission lines. 



0602+673 



1438+385 
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Observations at the Keck I Telescope by Taylor show a spec- 
trum with broad Hell at 4538 (z=1.766), broad Cm] at 
5311 (z=1.783), weak Aim at 5133 (z= 1.762), and weak Mg ii 
at 7783 (z=1.780). We estimate the redshift to be 1.775+0.01 
based on the stronger lines. Earlier, Vermeulen et al. (1996) had 
published a probable emission line redshift of z= 1.773. 

1459+480 

The emission line redshift, z=l. 059+0. 002, is based on broad 
Mg 11/12798 and narrow [Oll]/l3727 emission lines in an un- 
published Keck I Telescope observation on 1996 June 18 by 
Vermeulen & Taylor. 

1716+686 

The probable redshift, z=0. 339+0.001, is based on an observa- 
tion at the Keck I Telescope on by Taylor and strong [O ll], H/?, 
[Olll], and Ha. This redshift is inconsistent with the z=0.777 
found by Kuhr (1980) and also with the z=0.798 published by 
Stickel&Kuhr(1994). 

1734+508 

The emission line redshift, z=0. 835+0. 001, is based on very 
prominent narrow [Nev]/13426, [Oll]/i3727, [Nelll]3869, and 
[O lll]/14959, 5007 emission lines in an unpublished Keck I 
Telescope observation on 1996 June 18 by Vermeulen & Taylor 

1809+568 

The emission line redshift, z=2. 041+0.002, is based on promi- 
nent ClV/tl549 and C III] /i 1909 emission lines in an unpub- 
lished Keck I Telescope observation on 1996 June 18 by 
Vermeulen & Taylor. 

2054+611 

The emission line redshift, z=0. 864+0. 002, is based on promi- 
nent Mg 11/12798, Hy/14340, and J^/14861 emission Unes in an 
unpublished Keck I Telescope observation on 1996 June 19 by 
Vermeulen & Taylor [O Ill]/14959 is probably present at the edge 
of the observed spectrum as well. 

2229+695 

The emission line redshift, z=l. 413+0. 002, is based on broad 
Cm] /1 1909 and Mgll/12798 emission lines in an unpublished 
Keck I Telescope observation on 1996 June 19 by Vermeulen & 
Taylor 

2238+410 

The emission line redshift, z=0. 726+0. 001, is based on narrow 
[Oll]/i3727 and [Olll]/14959, 5007 emission lines in an un- 
published Keck I Telescope observation on 1996 June 18 by 
Vermeulen & Taylor. Ca H+K absorption at the same redshift 
is also visible. 

2319+444 

The redshift, z=1.251+0.003, is based on a broad Mg 11/12798 
emission line and a clear 4000 break in the continuum, seen in 
an unpublished Keck I Telescope observation on 1996 June 19 
by Vermeulen & Taylor. 



Appendix B: Kinematic modeiing procedures 

The starting point in all cases was the set of circular Gaussian 
parameters estimated from the complex visibilities per epoch as 
discussed in Paper I. These include the positions {x,y) of the 
components found by difmap with respect to the reference point, 
their uncertainties (cr^, a-y), and the elements of the correlation 
matrix associated with these parameters. In the ideal case, we 
could estimate the kinematic parameters for all of a source's 
components in a single fit, incorporating the difinap-ouVpul cor- 
relation matrices as a priori covariance matrices to reflect the 
fact that the visibilities and their (m-v) distribution may constrain 
different parameters more or less effectively in different epochs. 
This approach should preserve most of the "information" con- 
tained (or not) in the original visibilities. 

However, the existence of a "splitter" component, where two 
distinct difinap components are identified with a single jet com- 
ponent at some epoch(s), can frustrate this approach. The resid- 
uals of the individual "sub-components" with respect to the es- 
timated position for the associated single component dominate 
the chi-square statistic for the kinematic model, which in turn 
becomes insensitive to adjustments in the parameters of all other 
components. Additionally, the use of an a priori covariance ma- 
trix that includes non-zero cross-component elements ensures 
some of the specious adjustment to "splitter" components com- 
municates to the other better-behaved components. Thus, the 
first variant approach to the kinematic modeling was to fit inde- 
pendent models for each component within a source rather than 
for each source as a whole. This meant not taking into account 
any cross-component correlation-matrix elements in forming the 
a priori covariance matrices. Further, the effective position of 
a "splitter" at an epoch was controlled by the relative sizes of 
the constituent sub-components' uncertainties. In some cases, 
the resulting time series of component position was not well- 
described by monotonic motion. We also tried other weighting 
schemes for the (ctv, cr,,). Each of the above variants became an 
independent choice in the tactics for estimating the kinematic 
models for the sources: 

- include all components in a source in a single fit or fit inde- 
pendently per component 

- use the difinap correlation matrix as an a priori covariance 
matrix or not 

- incorporate (ctv, cr,.) as estimated from difmap, treat them as 
uniform across all epochs, or set them proportional to recip- 
rocal flux-density 

The complete set of kinematic models could be computed for all 
sources given a combination of the above three choices. 

In order to gain confidence in the incorporation of "splitters" 
into the statistics of component kinematics, we had to investigate 
whether using any of the alternate tactics led to any biases. First, 
we looked at only sources composed of Ql components. These 
could be fit using the ideal method as described above. We then 
recomputed the kinematics fits for these Ql-only sources using 
different tactics more accommodating to "splitters". Specifically, 
we would thus be comparing separate fits using the first and the 
last option in the each of the above three bullets. Let us call 
these the "Joint" and "Independent" fits. To be able to include 
the lower-quality components reliably, we would need to demon- 
strate at least that the "Independent" fit introduces no bias com- 
pared to the more rigorous "Joint" fit for these Ql-only sources. 

Fig. IB. II shows the correlation between the proper-motion 
estimates for these two sorts of fits - panel (a) for fi" and panel 
(b) for yu''. Each point corresponds to a single component, whose 
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"Joint" vs. "Independent" Fits 
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Joint Fit [mas/yr] 

(a) 

jJ: "Joint" vs. "Independent" Fits 
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Fig. B.l. Comparison of proper motion values computed using 
different tactics in the kinematic modeling. Each point represents 
a single component, with uncertainties along the abscissa from 
the "Joint" fit and along the ordinate from the "Independent" fit. 
Separate 1-D absolute-deviation fits are shown, one using the 
"Independent" fit as the dependent variable (red) and the other 
using the "Joint" fit as the dependent variable (blue). The dotted 
green line marks the (ideal) diagonal. Panel (a) shows and 
panel (b) shows /v''. 



uncertainty in the abscissa direction come from the "Joint" fit 
and in the ordinate direction from the "Independent" fit. We 
made two separate one-dimensional absolute-deviation straight- 
line fits to these points, one using the "Independent" fit as the 
dependent variable (red) and the other using the "Joint" fit as 
the dependent variable (blue). Fit statistics are annotated in the 
appropriate color. A dotted green line marks the diagonal, upon 
which all points would ideally lie. We then computed the per- 
pendicular distance on these plots from each point to the diago- 
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Fig. B.2. Histograms of yUj^ normalized by o"^^ . Here, jt/^ is the 
perpendicular distance from each point in the panels of fig- 
ure lB.ll to the diagonal. The top panel shows ju^ and the bottom 

shows yUj^. 

nal (lyUjoint - yUindepI/ V2) and the uncertainty in this jt/j_ via con- 
ventional error propagation assuming no correlation between the 
"Joint" and "Independent" uncertainties. This perpendicular dis- 
tance is significantly more strongly peaked near zero than is ei- 
ther of the uncertainties in the two kinds of fits, and there are 
no obvious trends in the two-dimensional plot of fi^ vs. either 
uncertainty (here with keeping a sense of above or below the 
diagonal). Fig. IB .21 shows histograms of //^ normalized by cr^^ 
- distributions that are consistent with 0-mean processes, signif- 
icantly more compact and peaked (kurtosis = 21.2 for 11.8 
for //^) than Gaussian. 

Thus, we should be able to use the "Independent" fit tactics 
in order to be able to incorporate components from sources hav- 
ing lower-quality ones without biasing the statistical results from 
the higher-quality sources. We scale the resulting uncertainties 
from these fits to a reduced-chi-square of unity, which should re- 
sult in reasonably conservative values. A disadvantage of this ap- 
proach is that some information contained in the difmap fits from 
paper I for sources containing only Ql components will not be 
carried forward through the kinematic modeling. Furthermore, 
the "Independent" fit tactics as described above produce the nu- 
merical result cr^v = cr^v in the course of the least-squares esti- 
mation, essentially circularizing the beam in all epochs. 

Appendix C: Kinematic model results 

Table IC.ll lists the results of the kinematic modeling for all 779 
components from sources participating in the proper-motion 
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analysis. Column (1) lists the lAU name, column (2) gives the 
optical classification (Q: quasar, B: BL Lac object, G: galaxy, 
U: unclassified object), column (3) lists the redshift, column (4) 
gives the jet component identification, column (5) the quality 
class as defined in section [3] The following three columns (6), 
(7), and (8) give the number of times this component has been 
detected, the time span from the first to the last detection, and 
the reference epoch. Column (9) and (10) list the component's 
position at the reference epoch in rectangular coordinates. 
Columns (11) and (12) list the proper motion components in x 
and y respectively, and the next pair of columns (13) and (14) 
show the radial and orthogonal proper motion components. 
Column (15) gives the total proper motion value, calculated 
from fi'\iJ.^'. The last three columns (16)-(18) list the radial, 
orthogonal, and total apparent velocities. Parameters in columns 
(9)-(18) are listed with their uncertainties. The uncertainties for 
the "splitters" are typically 2-3 times the uncertainties for a 
"normal" component, driven mostly by scaling the high reduced 
chi-square of their fit. Components that have observations 
at only two epochs do not have an associated uncertainty; 
with no degrees of freedom in the fit for the kinematic model 
(four constraints and four unknowns), the scaling to x = 1 is 
ill-defined. In Table IC.ll we indicate this missing uncertainty by 
a bar. 

Additionally, we display the kinematic modeling results 
graphically . Fig. lC.ll shows the source 1 106H-380 as an example. 
We show three panels: 

Panel (a) shows a plot of the positions of all components at 
all epochs on a plane tangent to the core. The components are 
color-coded, and the epochs are represented by symbols (first 
= triangle, second = square, third = diamond, fourth = x, fifth 
= asterisk). The scaled uncertainties for each component from 
Paper I are also overplotted on top of the symbols. In addition 
to the observed positions plotted in color, the modeled positions 
at each epoch are shown in small black symbols. Each of these 
plots is drawn so that one mas on the sky has the same scale in 
both the X- and y-axes on the plot (the scale of course can vary 
from source to source). 

Panel (b) shows a plot of the x and y positions of the components 
as a function of time. The component color-coding and epoch 
symbols remain the same as in the X-Y plot. To guide the eye 
for small motions, a horizontal dotted line is plotted for each 
component. 

Panel (c) shows a plot of the total proper motion as well 
as the radial and orthogonal proper-motion components, for 
each jet component plotted as a function of radial distance 
from the core. The jet-component color-coding remains the 
same as in the other plots. The scaled errors in both r and the 
proper-motion component are plotted (unsealed errors used for 
components with only two epochs of observations). Plots for 
all 237 CJF sources are available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/AH-A/ 
In addition, the plots can be downloaded from 
a CJF-archive page at the MPIfR in Bonn 
( |http://www.mpifr- bonn.mpg.de/staf f/sbritzen/cjf.html[ ). 
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Fig. C.l. Example, for 1 106H-380, of the kinematic modeling figures. Panel (a) shows the positions of all components at all epochs. 
The components are color-coded and the epochs are represented by different plotting symbols. The modeled positions at each 
epoch are shown in small black symbols. Panel (b) plots x and y of each component as a function of time. For each component, a 
"stationary" dotted line is plotted for reference. Panel (c) plots fi'', and fi" as a function of radial distance from the core. The 
component color-coding of panel (a) is carried forward into the other two panels. 



Table C.l. The proper motion and ySapp-values of a subsample of the CJF sources (18 sources). The complete information for 237 sources is available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http:// cdsweb.u-strasbg.fr/c gi-bIri/qcat?J/A+A/| 
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